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- - - Abstract 

\o _ 

I We develop the symbolic representation method to derive the hierarchies of (2+ l)-dimensional integrable 

■ equations from the scalar Lax operators and to study their properties globally. The method applies to both 
(N ■ commutative and noncommutative cases in the sense that the dependent variable takes its values in C 

f~| . or a noncommutative associative algebra. We prove that these hierarchies are indeed quasi-local in the 

■ commutative case as conjectured by Mikhailov and Yamilov in 1998, |MY98) . We propose a ring extension 
based on the symbolic representation. As examples, we give noncommutative versions of KP, mKP and 

, Boussinesq equations. 
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1 Introduction 

Integrable (1 + l)-dimensional nonlinear evolution equations, i.e. equations of the form 

-^(^5 '^XX^ '^XX...X^ 

possess rich algebraic and geometrical structures such as Lax representations, a bi-Hamiltonian formulation, 
, ^ ■ infinitely many symmetries and conservation laws, etc. A lot of work has been devoted to the study of such 
equations and comprehensive classification results have been obtained. 



o 
o 



The extension of these remarkable structures to the (2 + 1) -dimensional case is not straightforward. It was 
proved in |ZK84j that there is no bi-Hamiltonian formulation for (2 + l)-dimensional integrable equations of 
, the same type as those for the (1 + l)-dimensional case like the Korteweg-de Vries (KdV) equation 

o • 

Ut = Uxxx + &UUx. 



In 1988, Fokas and Santini, cf. | SF88[ IFS88| . constructed a bi-Hamiltonian structure for the Kadomtsev- 
Petviashvili (KP) equation 

^ ■ Ut= UxxX + GUUx + ZD^^Uyy 

by considering it as a reduction of (3 + l)-dimensional system. Meanwhile, Magri and his coauthors, |MMT88] . 



> 

X 

■ explained this structure from the geometric point of the view by developing the concept of Nijenhuis G-manifolds, 
■ ■ ■ ' which amazingly works for both one and two space dimensions. Dorfman and her coauthors introduced the 
noncommutative ring of formal pseudo-differential operators, cf. (DF921 rA.D93| . They proved that the Fokas- 
Santini bi-Hamiltonian structure of the KP can be obtained from the Adler-Gel'fand-Dikii (AGD) scheme by 
considering the second order Lax operator 

L = Dl+u~ dy. 

The bi-Hamiltonian structure naturally leads to a recursion operator for the equation. However, the operator 
(so-called operand) takes its value in the ring of formal differential operators and is of a different type from the 
one for KdV (and therefore does not contradict the result of |ZK84p . A hierarchy of infinitely many symmetries 
can only be produced by making suitable combinations of the operator acting on distinct seeds jMMT88] . 

One of the main obstacles to extend the spectacular results of (1 + l)-diniensional integrable equations to the 
(2 -|- l)-dimensional case is that the equations themselves, their higher symmetries and conservation laws are 
non-local, i.e. the appearance of the formal integral or Dy"^ . In 1998, Mikhailov and Yamilov, |M Y98j . 
introduced the concept of quasi- local functions based on the observation that the operators and 
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never appear alone but always in pairs like D^^Dy and Dy^D^ for all known integrable equations and their 
hierarchies of symmetries and conservation laws, which enables them to extend the symmetry approach of 
testing integrability |MSS91j to the (2 + l)-dimensional case. 

In this paper, we develop the symbolic representation method to produce hierarchies of (2 + l)-dimensional 
integrable equations from scalar Lax operators and to study their properties globally. The method applies 
both to the commutative and noncommutative case in the sense that the dependent variable takes its values 
in C or in a noncommutative associative algebra. We prove that these hierarchies are indeed quasi-local in the 
commutative case as conjectured by Mikhailov and Yamilov in 1998, ^ .MY98, , . This concept of quasi-locality 
has to be extended in the noncommutative case. Here we propose a ring extension based on the symbolic 
representation. As examples, we give noncommutative versions of the KP, mKP and Boussinesq equations. 



2 Quasi-local polynomials and symbolic representation 

The symbolic method was introduced by Gel'fand and Dikii in 1975 |GD75| . However, its basic idea can be dated 
back to the middle of 19th century. Recently, we successfully applied this method to the classification of both 
commutative and noncommutative (1 + l)-dimensional homogeneous evolution equations in a series of papers, 
cf. |BW98, BSW98, OWOO . With the help of number theory, it enables us to give a global description of their 
integrable hierarchies |Wan 98 . The symbolic method is also powerful in dealing with differential (cf. |SW02| ') 
and pseudo-differential (cf. [MN02j ) operators. The authors of |MN02j generalised the standard symmetry 
approach MSS91 and made it suitable for the study of nonlocal and non-evolutionary equations, see also 

|MNW06| .^ 

In this section, we'll extend the symbolic method to the case of two spatial variables x and y. For simplicity, 
we restrict our attention to a single dependent variable u. Extensions to several independent variables and 
dependent variables are straightforward. 



2.1 Quasi-local polynomials 

We begin with basic definitions and notations of the ring of commutative and non-commutative differential 
polynomials. The derivatives of dependent variable u with respect to its independent variables x and y are 
denoted by Uij = dl.d^yU. For smaller i and j, we sometimes write the indices out explicitly, that is Uxxy and u 
instead of U21 and woo- A differential monomial takes the form 

We call k the degree of the monomial. We let denote the set of differential polynomials of degree k. The 
ring of differential polynomials is denoted hy U = ®k>il^^ ■ Notice that 1 ^ U and it is a differential ring with 
total x-derivation and y-derivation 

-t-OC -t-OC +OC +00 



Dx^^^u.+i^,— and Dy = ^Yl 



i=0 j=Q '■J i=0 j=0 'J 

The ring U is commutative if the dependent variable u takes its values in a commutative algebra, for example 
the ring of smooth hmctions in x and y. Let us denote 

e = D-'Dy, e-' = Dy'Dx. (1) 

The concept of quasi-local (commutative) polynomials U {Q) was introduced in |MY98| to test integrability 
of a given equation. To define it, we consider a sequence of extensions of U. Let QU = {Qf : f G U} and 
Q-^U = {0-1/ : f eU}. We define UoiQ) = U and Uk{e) is the ring closure of the union 

Uk-i{e) u euk-i{e) u e-'Ut^e). 

Here the index k indicates the maximal depth of nesting the operator in the expression. Clearly, we have 
Uk-i{&) C Uk{Q)- We now define U{Q) = limfe^oo ^*:(9)- However, for a given / e U{Q), there exists k such 
that / G Uk{&)- The extension U{Q) has a natural gradation according to the number of u and its derivatives 

Note that lAk{&) is not invariant under transformation of variables. For example, a simple transformation 

x' ^x + y, y' ^y (2) 



2 



transforms 9 = D^^Dy, ^l-Q = 1- D-^Dy. Hence 9"^ (1 - 9)"\ which is not in Uk{Q). 

If u takes its value in a noncommutative associative algebra, the ring lA is noncommutative. Typical examples 
are the algebras of n x n matrices and Clifford algebras (see |OS98| for more examples and noncommutative 
(l + l)-dimensional integrable evolution equations). For any element /, g and ft, in a noncommutative associative 
algebra, we use the notation 

Lf{g) = fg, Rf{9)^gf (3) 

for the operators of left and right multiplication and 

Cf{g)^Lfig)-Rf{g) = fg-gf (4) 

for the commutator (notation adf is also commonly used). It is a derivation since it satisfies the Leibniz rule: 
Cfigh) = Cfig)h + gCf{h). 

So far there is no proper extension from noncommutative differential polynomials to quasi-local noncommutative 
polynomials, mainly because the study of noncommutative (2 + l)-dimensional integrable equations is a new 
and challenging topic. We can take the corresponding extension of Z//(9) as the starting point. However, this 
turns out to be too restrictive. This will be discussed further in Section 



2.2 Symbolic representation 

The basic idea of the symbolic representation is to replace Uij by u^^^rj^ , where ^ and are symbols. Now the 
total differentiation with respect to x, that is, mapping utj to Ui+ij, is replaced by multiplication with ^, as in 
the Fourier transform. Similarly, the total differentiation with respect to y mapping Uij to u^j+i, is replaced 
by multiplication with 77. For higher degree terms with multiple it's, one uses different symbols to denote each 
of them. Its symbolic representation depends on where u takes its value, i.e. whether it is commutative or not. 
For example, the noncommutative binomial UijUki has symbolic representation u'^£,iiliS,2''l2- the commutative 
case, one has UijUki = uuUij. We therefore need to average its symbolic representation over the permutation 
group E2 so that UijUki and UkiUij have the same symbolic form. So the symbolic representation of UijUki is 
■t(^i'^i^2^2 + ^IviiiVi)- Symmetrisation makes the symbol representation of monomials unique. 
Let A'^ be the space of algebraic polynomials / in 2k variables, and rji, where i — 1, ■ ■ ■ ,k. 

The symbolic representation defines a linear isomorphism between the space U'^ of (non)-commutative differential 
polynomials of degree k and the space A'' . It is uniquely defined by its action on monomials. 

Definition 1. The symbolic representation of a differential monomial is defined as 

I u'' Ci rji^ £,2^ vi^ ■ ■ ■ Ck Vk" ^ (noncommutative); 

Ti^l . . . f^k n^l' 
^cr(l)'/(T(l) ^a(k)'l<y(k)^ 



1 TrE.e..Ci)^^a)---C.)<.)- {commutat^ve). 



There are two parts in the symbolic representation of a monomial: the first part, , indicating its degree; the 
second part being in A^ . The symbols of Gel'fand and Dikii only contain the second part. We emphasise that 
the first part is very important in dealing with the monomials the case of several noncommutative dependent 
variables. 

In general, we denote the symbolic representation of P S ^ whether it is commutative or not, by P and 
Q hy Q. The multiplication of P and Q corresponds to the symbol 

PQ{£i,mr-- ,£k,vk,£k+i,vk+i, - ■ ■ ,£k+i,vk+i) = P{£i,rii,--- ,£k,Vk)Qi£k+i,Vk+i, - ■ ■ ,£k+i,vk+i); (5) 

when P and Q are commutative differential polynomials, the right-hand side needs to be symmetrised 
PQiCvi,- ■ ■ ,£.k,ilk,£k+i,Vk+i, - ■ ■ ,£k+l,Vk+l) 

^ (j^ _)_ /)! X/ Pi^a{l)^Va(l),- ■ ■ ,£a(k),Va(k))Q{£cr(k+l),Va{k+l)i- ■ ■ ,£a{k+l),'na{k+l))- (6) 

The total derivatives and Dy have the following representations: 

DxP(.£.i,m,£2,m, ■■■ ,£k,vk) = {£i-\ hCfc)P(Ci,m,6,'72, • • • ,£k,vk); 

DyPi£i,m,£,2,m, ■ ■ ■ ,£k,vk) ^ {m + — \'Vk)P{£,i,vi,£.2,m, ■ ■ ■ ,£k,vk)- 
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Naturally, the actions of O = ^Dy and 8 ^ on P E U'' can be represented as 

ryi H Vrjk 

^^^Pi^i,Vi,^2,V2, ■ ■ ■ ,£.k,ilk) = ^-T^PiCi,Vi,£.2,V2, ■ ■ ■ ,ik,ilk). 

t.1 + ■ ■ ■ + t.k 

Similar to the multiplication of P and Q, we have 
Definition 2. Let P e U'' andQ eUK Then 

PBQ = ■ • • ,6,??fc)^^^-^-4 "t Q{^k+i,Vk+i, ■ ■ ■ ,^k+i,Vk+i)- 

H 1- Vk+l 

When P and Q are commutative differential polynomials, the right-hand side needs to be symmetrised. 

In the same way, we can uniquely define the symbol of PQ~^Q. Together with the multiplication rule lO or 
we have now completely determined the symbolic representations of the elements in Ui (O). By induction, we can 
define the symbolic representation of any element in U{0), which is a rational function with its denominator 
being the products of the linear factors. The expression of the denominator uniquely determines how the 
operator 8*^ is nested. 

Example 1. The expressions u^^i '^^^^'^ ~| (^nd '^^Ci^^ '""g symbolic representations of UxQ{uQ^^u) G 
^^2(8) and Ux{Qu)Q^^u Ui{Q) respectively. 

Let us define the symbolic representations of pseudo-differential operators motivated by their Fourier transforms. 
First we assign a special symbol X to the operator D^- Its formal inverse has the symbol Then we 
have the following rules for / e U^{Q) with symbol / and g £ U\Q) with symbol g: 

X o /(a, '71, . . . , 6, ^fc) = (Ci + • • • + ik + • ■ • , ^k,Tik), 

which represents the Leibniz rule ° f — Dx{f) + / ° Dx] and the composition rule 

fXP o gX^ ^ /(Ci, 771, . . . , ^k,VkMk+i + ■■■+ (k+i + Xfg{^k+i , Vk+u • • • , ^k+um+i)X''- 

This composition rule is valid for both positive and negative powers p and q. For positive powers, it is a 
polynomial of X. For negative powers, one can expand it at X = 00 to identify it. 

Example 2. The symbol of D^^u is if we assign ^ as the symbol for u. Then in the noncommutative case 
we have 

D^^uD^^u = u^D^^ — {2uux + Uxu)D~^ + {"iuuxx + + Uxxu)D7^^ + 

Symbolically, we can compute it by 

u u 2^ 1 6+26 , e? + 366 + 3e| , , 

o = -, r = U ' ' ^ I \ 



x + i x + i + + x^ 

To extend the symbolic representation from one dependent variable to several noncommutative dependent 
variables is straightforward. We need to assign new symbols for each of them such as assigning rp-^ for u 
and ^(^^jTy^^) for v and so on. As one can imagine, the abstract formulas for arbitrary n dependent variables 
get very long and it is difficult to present them in a compact way. We have no problems in handling small n or 
given expressions. 

Example 3. The symbolic representation of UxVy is uv^'^^^ rj'^'^^ and that of VyUx is vuS}^^^ ri'^'^\ 

Notice that the only difference in the symbolic representations for UxVy and VyUx is uv instead of vu, that is, 
the non-commutativity of dependent variables lies in the first part of the symbolic representation indicating 
the degree of the expression. 



3 Lax formulation of (2 + 1)— dimensional integrable equations 

In this section, we give a short description of construction of (2+l)-dimensional integrable equations from a given 
scalar Lax operator based on the well-known Sato approach. For details, see the recent books jBSOll |KupOO| 
and related references in them. 
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Consider an m-th order formal pseudo-differential operator 

A = a^D™ + ara-iD'^^^ H h ao + a^iD^^ H , m > 0, 

where the coefficients are functions of x and I?^^ is the formal inverse of the total x-derivative D^. Here 
the functions take their values in commutative (for instance complex numbers C) or noncommutative (for 
instance the algebra of rt x n matrices) associative algebras. 

To define the product of two pseudo-differential operators requires the action of differential operator £)" on a 
multiplication operator / (given by a function of x): 

Let the commutator be the bracket on the set of pseudo-differential operators. Thus, the set of pseudo-differential 
operators forms a Lie algebra. For an integer < to, we introduce notations: 

A>k = flm^T + a„_iL>™-i + • • • + akDl 
A^k^A- A>k - ak-iD';-^ + ■■■ 

When k G {0, 1}, the Lie algebra decomposes as a direct sum of two subalgebras in both commutative and 
noncommutative cases. We denote the projections onto these subalgebras by V±. Such decompositions are 
naturally related with integrability and lead to admissible scalar Lax operators for (1-1- l)-dimensional Lax 
dynamics (cf. |Bla98j pp. 270 for commutative case), namely 

k = 0: L = -Fu("-2)£)j-2 _^y(n-3)£,n-3 ^ ... ^y(o)^ ^ > 2; 

k = l: L = + u^''~^'>D^~^ +u^''-^'>D^~^ + ■■■ + u'^'^'> + D-^u'^-^\ n > 2; 

where u^*' are functions of a spatial variable x. 

We know every pseudo-differential operator of order n > has an n-th root. This allows us to define the 
fractional powers L^,i S N. Let Bi = 'P+(L"). For each choice of z, we introduce a different time variable U. 
Thus, the fiows defined by 

df 

^^ = [B.L], .en 

commute, 

This setting up can be generalised into the (2-|- l)-dimensional case as follows: An m-th order pseudo-differential 
operator of two spatial variables is of the form 

H = -Dy + amD"^ + ara-iD'^-'^ H VaQ + a^iD^^ H , to > 0, 

where coefficients Ofc are functions of x, y and for an integer k < m, we split into 

H>k = ttmD"^ + a„i_il?™ ^ + • • • + cikD^ 
H<k = H- H>k = -Dy + ak-iD^,~^ + ■■■ 

Similar to (1 -I- l)-dimcnsional case, this operator algebra decomposes as a direct sum of two subalgebras in 
both commutative and noncommutative cases when k G {0, 1}. The admissible types of scalar Lax operators 
for the case of (2 + 1) dimensions are L — Dy, i.e. 

a. fc = 0: L = £)«-2 ^ y (n-3) £,«-3 _^ y (0) n > 2; 

b. fc=l: L = D^' -f -hu("-2)£)^'-2 + . .. + ^(0) „>2; 

c. fc=l: L = u(0) -l-D-iu(-i) -Dy; 

where u*-'^ are functions of two spatial variables x, y. We often use u,v,w, - ■ ■ in the examples. When n = 1 
for case b, we have L = D^ + u*^"^ + D^^u^~^^ — Dy, which can be transformed into case c by transformation 
(0). Therefore, we exclude it from our consideration. 

Let S = Dx + ao + a^iD^^ For any operator L as listed in cases a, b and c, the relation 

[S, L] := SL-LS^ 0, (7) 
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uniquely determines the operator S by taking the integration constants as zero. Furthermore, we have [S*", L] = 
for any ri G N. For each choice of i, we introduce a different time variable U and define the Lax equation by 

I^I^V.i.], (8) 
where k is determined by the operator L as listed in cases a, b and c. 

Theorem 1. For the operator S uniquely determined as above by the flows defined by Lax equations ^ 
commute, that is, dtjOt^L = dt^dt^L. 

Proof. Using Lax equations, we have 

dt^duL = [dt,S'>u.L\ + [S\k,dt^L] = [dt^S^u.L] + [5%,., [5^>fc,L]] 

Hence 

Meanwhile, from formula (|7|), we get 

[duS, L] = ^[S, duL] = -[5, [S\k. L\] = [[5Vfc, S],L\ 



implying 



Thus 



duS^[S\u.S]. (9) 



dt,S'>k-duS-'>k-[S'>k.S\k] 

= — [<S'-'<fc, 5*>fe + 5*<fe]>fe — [S^>k, S^>k + S'-'<fc]>fc — [S-' >k, S^>k] = — ['S'-'<fe, 5*<fc]>/c 
= 0. 

This leads to dtjdt^L = dtidt^L. o 

Remark 1. In the commutative case, there is one more admissible type of scalar Lax operator, that is, 
k = 2: L^ u^'^^^D: + u("-i)i?r ' + u'"-')i?r' + • • • + + D-\^-'^ + D^^u^-^^ - Dy. 
Theorem^ is also valid for this type by taking S — u^^^^ + oq + a_i_D^^ + • ■ • instead. 

Example 4. Consider the Lax operator L — u^D^ + vD^ + w + D^^p + D^^q — Dy, where all dependent 
variables take their values in C. Using formula Qj, we have S — uD^ + ag + a-iD^^ + ■ ■ ■ with 

ao = --;^(Ux h B-) 

2 u u 

a_i = -i2u,, - ^ — - + J - ^ + -D-^Dy— + e -e- + - e-f . 



and S^>2 = 0, S^>2 = u^D^, S^>2 = u^D^ + 3u^(mj; + ao)D^. By the Lax equation, we obt 



am 



""ts = jiu^Uxxx " Su^ei - 3u{Q^)y + 3u^Wj;a; + Su'^Wx + ivy + 3vVx - 3^1^^ - 6{uv)xQ^) 

ft3 = U^Vxxx + Su'^iWxx + UPx + U^p + {v^x + 2Wx){u^ + fto)) 

Wt3 = u^Wxxx + iu^iwxx +Px){ux + ao) + 'ip{u'^aQ)x + "iu^iuqx - 2uxPx + "^Uxq) 
Pt3 = {u^p)xxx - 3 {[u'^{u.j,p + oqp + q)]xx + 2[u'^q{ux + ao)]x) 
qts = {u^q)xxx - 3[u^q(wx + aQ%.^ 

Notice that the reduction ofv~w=p = q = Q leads to well-known (2 + \)-Harry-Dym equation 

^t3 = -:{u^Uxxx - &UyQ- - 3m(6-)„) . 
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4 Lax formulation in symbolic representation 



In this section, we carry out the formalism of section in symboHc form. We first consider noncommutative 
dependent variables to simplify our formulae. The strategy for the commutative case is to do the calculation 
as much as possible by treating it as noncommutative and only do symmetrisation at the last stage to get the 
uniqueness of the symbolic representation since the symmetrisation complicates the calculation dramatically. 

Let us assign the symbols ^^^'^ry^*) for dependent variable m^*^ and the symbol Y for the operator Dy. The 
symbolic representations of the admissible scalar Lax operators are 

a. k = 0: L = X"-y + u("-2)j^n-2_^y(n-3)j^n-3^...^y(0)^ ^ > 2; 

b. fc = l: L = X"-y + M("-i)X"-i+u("-2)j^«-2^...^^(0)^^(-i)_i__^ ^>2; 



c. 



k=l: L = _r + u(")+u(-i) 



We first treat case a. It is convenient to consider formal series in the form 

j = 21=0 12=0 

where are functions of their specific arguments and jk is defined by the number of ik in the list of [ii, 12, ■ • • , i/], 
which implies that ji — 1 and jk > 1, fc = 1, 2, • • • , L For example, when ii — 12, the arguments of function 

It is easy to check that 

where the polynomial Ni is defined by 

I I 
Ni[x^,yuX2,y2. . . . , x,, y,; X) = + X)" - X" - ^ (12) 

i=l i=l 

Proposition 1. For any operator L in case a, if the formal series satisfies the relation [S, L] — (cf. 
we have for I > 1, 

a[""-''^=a,(^^-\,7(-\...,e]r\r;(-\x) = n(X+ ^ t^y^b,{^^\v^\ ■ ■ ■ ,^^\v^\x), (13) 

r=l s—r+1 

where the superindex is € {0, 1, 2, • • • , n — 2} and the suhindex jk is defined by the number of ik in the list of 
[ii,j2, • • • ,ik\- The function hi, I >1, is defined by 

with 

= &.-i(e\r;j-^...,cj:-\,j-\x+^(-))-6,.,(e^ ;>i (15) 

and the initial function ci{s}i\ri[ \ X) = 

Proof. We compute [L, S] symbolically and collect the coefficients of each algebraically independent monomial 
of the dependent variables. For the linear terms, we have for any i G {0, • • • , n — 2}, 

af i^' , ,X)^ Jf = (^« , ,X)X^. (16) 
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For the quadratic terms, when i ^ j, we need to compute 

= {X + Ci'^yx^ (61 , ryW , X + ) - 61 (Cp^ , , X)) 

and when i = j, we have 



This can be obtained by substituting ^[^^ = ^2*' and ri[^^ = vi^^ in A^^ . It follows that 

{ili2) I Ail) (il) Ai2) (12) Y\ — "^'^'^ 

2 ''Oi 'S2 '^/j2 '^'i /V J^i) fe) y^ 

= (X + e]r)r^X^^62(C]r\ry(:^\^2^\r7(:^\x). (17) 

Assume that formula (|13() is valid for degree I — 1. Let us compute a[*^*^ in the same way as computing the 
quadratic terms. Then we have 



p=2 r=2 s-r+1 

i-1 ; 



■r— 1 s— r+1 

; / 

r— 1 s— r+1 

; / 

nfy, \^ f(is)Yr^ /Ail) (il) Ail) (il) v^ 

r—1 s— r+1 

According to definition 114|l . we proved that formula is valid for degree I. o 
Similarly, we prove the following result for any operator L in cases b and c. 
Proposition 2. For any operator L in case h, consider a formal series of the form 

n— 1 n— 1 

S = X+Y^ ^.«a«(ef\r7«,X)+ ^^v^^4''^i^:\rl^:^,^\rl^,X) + ---,n>2 (18) 

2— — 1 ^1,^2 — ^1 

Here a.i are functions of their specific arguments, the superindex ZsG{ — 1,0, l,---,7i — 1} and the subindex jk 
is defined by the number of ik in the list 0/ [ii,i2, • • • ,?fc]. // it satisfies the relation [S, L] = (cf. we 
have for I > 1, 



(ili2---ii) _ ^ (Ail) (il) Ail) (il) y\_T1(u(a \, ST^ \V ^ /(^(^i ) A^i) Y\ (1Q\ 

r— 1 s— r+1 



where h{ir) ~ X + if ir = —1, otherwise h{ir) — X. The functions bi are defined in PropositionUl cf. 

formula and f J,?)) . 



8 



Notice that for operator L listed in case c, we have Ni{xi, yi,X2,y2, ■ ■ ■ ,xi,yi] X) = — Vii CDi which 
corresponds to the action of total y-derivative Dy. Thus 

Proposition 3. For operator L listed in case c, if a formal series of the form 

where the superindex is G { — 1,0} and the suhindex jk is defined by the number of ik in the list of [ii,i2, ■ ■ ■ , ik], 
satisfies the relation [S, L] — (of. we have for I > 1, the functions ai are defined in Proposition\^ cf. 

formula with Ni{xi,yi, X2,y2, ■ ■ ■ ,xi,yi]X) = - ^^Li Vi- 

To construct the hierarchy of the Lax equations we need to expand the coefScients of operator (|10|l or (|18|l or 
H2()|l at X = oo and truncate at the required degree. 

Definition 3. We say that a function h{xi,yi,X2, y2, ■ ■ ■ , xi,yi, X) is kth (quasi-local) polynomial if the first k 
coefficients of its expansion at X = oo are symbols of (quasi-local) polynomials. If a function is kth (quasi-local) 
polynomial for any k, we say it is (quasi-local) polynomial. 



When n > 2, that is, operator L in cases a and b, the expansion of Ni{xi, yi, X2, y2, ■ ■ ■ t^i, yi\ X) at X = cx) is 
of the form 

Ni{xi,yi,X2,y2, ■ ■ ■ ,xi,yi; Xy^ 



1^^^ iy( n \ y 



For operator L listed in case c, we know Ni{xi, yi, X2, 2/2, • • • , xi,yi;X) = — Vi- Therefore, if we want to 
prove that the coefficients of operators (|10|l . (|18|l and l|20() . i.e. the functions ai, are quasi- local, we need to 
show that the functions ci can be split into the sum of the image of and the image of Dy. It is clear from 
formula (|15|l that ai\S,i\ ri^\ X) are quasi-local since we have ci{^[^\ ri^\ X) = ^[^K Now we concentrate on 
the cases when ^ > 1. 

Proposition 4. The functions Q(^j*^\ 77^*^'', • • • , ^j*''', ryj*'"', X) for I > 1 vanish after substitution 

= f(^'-i) and n^''^ = -v^"''^ n^"~'^ 

?ji 'hi 'O2 '61-1 • 



We give its proof in the Appendix. In fact, this proposition does not lead to our intended conclusion that hi and 
thus ai are quasi-local since the objects are rational, not polynomial. For example, the expression — 
representing uQu — {Qu)u satisfies the above proposition. However, we can't write uQu — {Qu)u = Dxfi + Dyf2, 
where both /i and /2 are in U{Q). 

We first have a close look at the quadratic terms. For any operator L in cases a and b, let 77*^ = ^jr^jV Then 
the function cafC , r?!! , : , X), where ^Y, ^ and 7^ 0, is a polynomial of 9'^ , , r = 1,2. Notice that 

The affine variety (cf. |CL097| ) defined by -f and rf^^ + 77}^^ is 
From Proposition 21 it can be written as 

where both /i and /2 are polynomials of 6*^^,^]^, r = 1,2. From formula 121|l . the part that is not in the image 
of Dx or Dy only depends on 6*^^, r = 1, 2. 

For operator L listed in case c, since the function Ni{xi, yi, • • • ,xi, yi,X) = — Y^\^i yi, cf. Proposition we 
let = Sl''rf{ instead. Then the function C2(C!' , ^7 : CZ : v7 : where 77"^ 7^ and ry*^ 7^ 0, is a polynomial of 
djl,Cjl.^ f — lj2. The above discussion is vahd. 
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Setting ^^C^^O,we obtain the ftinction C2 {^f^ ,V^,Q,V%,X) equals 
1 1 



— for cases a and b; 91^ — for case c, (22) 

which is zero after we symmetrise it with the permutation group £2- This imphes that the functions 62 and 
thus a2 are quasi-local in the commutative case. 

4.1 The commutative case 

We first give the formulae to compute the high degree terms of operator 5* for the commutative case directly 
from the formulae we obtained in Proposition and |31 

We denote the list (ii,Z2,--- by / and Ii = {isi,is2r ' ' 1*5,-): where 1 < Si < S2 < ■ ■ ■ < Sr < I and 
1 < r < / — 1. We use I2 = {ipi, ip2:' ' ' j'i'Pi-r) denote the list by removing the elements of list Ii from list 
/. For any element v in the set {ii,i2,' ■ • we denote the number of v in the list / by =ffir- Then the 
coefficient of Yii g{i^i2 - ii}(^^''^'')'^''^ ™ operator s, cf. formula (|18|l and H2()|l . can be computed by 



n 

We denote it by 



(23) 



In particular, for any operator L in case a, we have 
^ ; I 



lljre{iii2---ii}^^ CTgS, r=l s=r+l 



For cases b and c, the expression {X + p) ^, where p — S,j + X]l=r+i appears in formula 119() . Its 
expansion at X = 00 is of the form 

-, +00 

whose coefficients are polynomials in p. This does not affect the proof and result. For case c, we need to 
exchange the symbols and 77^-^ in the proof as the discussion we did for quadratic terms. Therefore, we'll 
only give the proof for this case. 

The immediate consequence of the previous discussion on the quadratic terms is the following statement. 
Proposition 5. All quadratic terms in operator S are quasi-local. 

In fact, we can prove this statement is valid for alH > 0. 

Theorem 2. Every term in the operator S , cf. formula ^J^, and l^2()\) . obtained via Y2!^) is quasi-local. 



Proof. We proved the linear and quadratic terms in S are quasi-local, that is, the statement is true for 
Z = 1, 2. Assume it is also true for / — 1. By induction, the function dj^'^ cf. I|24|l . is polynomial of S^^J^^ and 

6*/^ = \ll ) iJT^i '' = lj2,---,Z — 1, which obviously satisfies Proposition 01 This implies that the function 

jni2---n jg ^Yie sum of the parts with factors <^j'^'' -I- ^j*^-* + • • • -I- S,^^^'^ or ^/j^ — Bj^ using the argument of the affine 
variety. From the recurrent relation in Proposition^ we can compute the factor for any — Bi^ (without 
symmetrisation), which equals 

f]{x + + . . . + d^-) + d'") + . . . + ^f'"- TT(x + + ■■■+ V''' 

- + +■■■+ ^r:' + 1:' + ■ • • + ^'r^)"- ri(^ + + ■ • • + ^'r^)^=' ■ 

q=l q=l 
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This expression is zero under the substitution S/, = f^'"^^ . . . ^ >■ — q and — E^^"^^ + • • • + f = 
implying we have either Ei^^Oj-^ — Oj^) or S/^ (0/^ — dj^) as a factor. We know 

(^/. - ^/j = ^/.(e^r^ + • • • + c^r^) - (^jr' + • • • + ^jr'); 

Thus the Z-th degree term of S is quasi-local, o 

This theorem implies that every terms in 5" is quasi-local. From Theorem ^ we have 

Theorem 3. The hierarchies of (2 + I )-integrable equation with scalar Lax operators are quasi-local. 

As mentioned in Remark^ there is one more type of scalar Lax operator in the commutative case. The same 
result can be proved by extending the symbolic approach used in [SWOO| for the (1 -f l)-dimensional case. 



4.2 The noncommutative case 



The extension of the concept of quasi- locality to the noncommutative case is rather complicated. Dx and Dy are 
the only derivations for the commutative differential ring. The extension simply enables us to apply D^^ and 
Z)~^ on the derivations. We know that the C/, cf. I@J, are also derivations for a noncommutative associative 
algebra and we need to take them into consideration. Based on Proposition ^ we propose the following 
extension using the symbolic representation. We start with ^ the space of polynomials in 2fc variables, and 
rji, i = 1, ■ ■ ■ , k. Let — Xzi and rji = ^Zi, where A, /i are constants. We define I{A'^) as a set of all / € A'' 
satisfying /^_|_ ^ _q — 0. For the elements in 5 e I{A''), we compute ^^^.^.^^^ and and collect them 

in the set £xt{A''). Now we define a sequence of extensions of ^ = UkA'^. Let ylg^* — A and Af^* be the 
ring closure of the union Aq^* U {Lik£xt{A'')). The number of variables is a natural grade on ^f^* and again 
denote AS^*^ the space of elements in A^^ with 2fc variables, and rji, i — 1, ■ ■ ■ , fc. In general, we can define 
Af^* is the ring closure of the union AlZl U {Uk£xt{Af^\ )). Clearly, we have ylff* C ylf"^*. We now define 
A"""* = lim/^oo^r*- 

The above ring extension is bigger than we require. For example, the expression ~ is not equal to zero 
under the substitution ^1 -|- ^2 + = and rji + r]2 + rjs = 0, that is, not satisfying Proposition ^ But it 
is in I[Af^* ). If we start from the symmetric polynomial, we hope to end up the extension of commutative 
differential polynomial, i.e. quasi-local polynomials we defined in section[21 It is not difficult to notice that this 
extension is bigger than the concept of quasi-local commutative polynomials. However, neither can we formulate 
the such extension using symbolic representation nor starting from noncommutative differential polynomials as 
in commutative case. It is still an open problem. 

Here we compute explicitly some Lax flows for lower order scalar Lax operators of noncommutative dependent 
variable, which will help the reader to see the structures. 



4.2.1 Noncommutative Boussinesq equation 



Let us compute the hierarchy of Lax operator L = + uD^ + v — Dy. From Proposition ^ we have 



{X + - X3 - ryf ) {X + e("^)3 - X3 - ryf ) 

U\X + S}^^)X 



+ 



+ 



+ 



{X- 




-|-e«)3-X3- 

uv{X + ^^°'^) 




- 4" 








)3 - (X - 






(X- 






X3 




{X- 




4.^(0))3_X3- 

vuX 


1 ^\ 








nr. 


)3 - (X - 

AO) 






(X- 


Hi"' 




X3 





iX ^ 3X2 + ( ^ + ^ T^X3 
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3 2 

-uC['^ 2u7^['^ 2<<"^ t^i; w aje^'^ , 5$^^ , 1 

^^9 9^9 ^9^(0) 9 9 ^"^27 ^ 27 ^ ^1'^+^^'' JJ^4+---' 

which corresponds to 

2uy 2vxx ©i" i^'^' "^^u 7uUa; 5uxU D^^{{Qu)u) D^^{uQu) _^ 
9 9~'^~9~^~9 9 9"^^T~^^T~^ 27 27 ' ^ 

This leads to 

Sio = Dt + ^uDl + ^(u, + 2v)Dx + ^Qu + ^u^ + ^u,, + ^w,; 

5 5 5 5 10 5 

S>o ^Dl + -uDl + -{ux+ v)Dl + {-Qu + -u" + —u^x + ■^Vx)Dx 

10 5 5 

+ -^Vxx + ^(©'y + uv + vu) + — {uux — UxU + D^^{{Qu)u — uQuj). 

And we have 

Ut2 = -Uxx + 2Vx 

Vt2 = Vxx + l{Uy - U3x - UUx +UV - vu) 

EHminating the dependent variable v from this equation and writing t2 — t, we can derive the (2 + l)-dimensional 



1 


4 


2. 2 


2 




Utt = -^U4x H 






)xx H~ "^^x 


[u, D^^ut]. 


transformation x ^/6x, 




' fy, u 1 




■ 2t, it can be rewritten as 


Utt = -Uix 






CX ~l" ["^7 


Dx'ut]. 



Let w — ^u. We get (2 + l)-dimensional noncommutative potential Boussinesq equation 

Wtt = -W4x + Wxy - Hwl)x + [Wx, Wt] 

4.2.2 Noncommutative KP equation 

Let us compute the hierarchy of the Lax operator L = + u — Dy. Since the order of L is low, it is easier to 
compute the operator S directly order by order as in Exercise 01 than to apply Proposition |2 Let 

S^Dx + a^iD-^ + a^2D-'^ + a^^D-^ ■■■ . 

Using formula pO(l . we obtain 



This leads to 







Q-l 


^ 2 








a-2 


_ Ux 

4 








a-3 




2 






8 


~8~ ^ T ~ "s" 


'^>0 




3u ^ 
^^Dx 


+ -^[Ux 




•^>0 


= Dt- 


h 2uDl - 


f {2ux 4 





"" 2 2 
3u G^u 1 

2DlDy -Dl + 2uDy + QuDx + ^ + — - l^D-^CuQu + L^. 
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And we have 

Ut3 = ^{Uxxx + iuux + iUxU + iQuy - 3C„9u), 

which is the noncommutative KP equation and can be obtained from the recursion operator of KP, of. |DF92 | , 
and 

""44 — lj^{''^xxy + Q^Uy + 2uUy + 2UyU + UxQu + {Qu)ux + CuD^^Cu&U — QCuOu — CuQ^u), 

which appeared in jHTn8| . cf. formula (2.14), with a mistake: the last two terms are missing. 
4.2.3 Noncommutative mKP equation 

Let us compute the hierarchy of the Lax operator L = + uDx + v + D^^w — Dy. Let 

S ^ Dx + ao + a^iD-^ ^ . 



Using formula pOfl . we obtain 



This leads to 

S>i^Dx; Sl^ = Dl + uDx\ Sl^ = Dl + + -{2ux + + 4v + 2{Dx + ^)-^Uy)Dx. 



And we have 



Ut2 = Uy + 2Vx + uv — vu 
Vt2 — Vx + 2Wx + UVx + uw — wu 
Wt2 = -Wxx + {wu)x 

-\C.^Ux + iwx + f C„w + KZJy ~Dl- uDx + Ru^ - a,)(Dx + ^V^Uy 

Vta = Vxxx + ^UVxx + jUxVx + fw^Ul + fwi + ^UWx + f Wx" + jWUx + jUxW - fCtyU^ 

~2CwV + l{Rv^ ~ Cyj){Dx H ^Uy 

^ Wts = Wxxx - {wao)xx + (wa-i)i 

The reduction v = w = leads to the noncommutative mKP equation 

1 3 3 3 3 3 C 

= -Uxxx + -^CuUxx - -^uuxu + -Uxy + + -(^^y " -D^ " LuDx + RuS){Dx + -;7-)~^Wy (26) 

If we introduce new variable p satisfying Uy — Px + ^CuP, The system for u and p is equivalent to the matrix 
mKP in IXdil. 



If u does not depend on y, i.e. Uj, = 0, this gives us the noncommutative mKdV equation, j()S98l [r)WOO| 

1 3 _3 

'^is — ."^xxx ~^ Q^u^xx Q^^x'^- 
4 o o 

If u takes its value over C (commutative), this gives us the well known mKP equation 

13 3 3 

Ut3 = -Uxxx - -^u'^Ux + -<duy + -UxQu. (27) 
4 o 4 4 

We know Miura transformation V = Ux — \ + Qu transforms the mKP equation (|27(l into the KP equation 
Vt3 = jVxxx + jVVx + |6Vy. However, we do not know a Miura transformation for noncommutative mKP 
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4.2.4 Nonsymmetric Novikov-Veselov equation 

Consider the Lax operator L = u + D~^v — Dy (case c) ^. Let S = + uq + a^iD^^ + • • • . Using formula 
(0, we have 

a_i = {Cu - Dyy^{vx + C^ao), 

leading to S^>i = D^, S^>i ^ Dl + 2aoDx, S^>i ^ + ia^Dl + 3(aox + a-i + al)Dx. From the Lax 
equation, we obtain 

f = + 2vx + 2aoUx + 2aow - 2uao 

+ 3{aQUx)x + 'ia-iUx + ialux + SaoVx + 3(wao)x — 3Cv{aQx + fl-i + a§) 
\ = - 3{vxao)x + 3{va^i)x + 3(val)x 

The reduction v ^ leads to 

Ut2 = ""xa; - SKDj^ - Cuy^UxjUx- 

Let w — {Dy—Cu)^^Ux. This equation transforms into the noncommutative Burgers equation — Wxx—2wWx, 
which is linearisable, that is, we obtain pt = Pxx by pw = —px- 

There is no reduction u = in noncommutative case. In the commutative case, the reduction u — leads to 
the nonsymmetric Novikov-Veselov equation jNV86| 

Vts = Vxxx - 3{vQ'^v)x. 

5 Conclusion 

We have demonstrated the power of the symbolic representation in the study of (2 + l)-dimensional integrable 
equations. It enables us to produce the hierarchies of (2 + l)-dimensional integrable equations from the scalar 
Lax operators and to study their properties globally. We proved the conjecture of Mikhailov and Yamilov on 
the ring extension for (2 + l)-dimcnsional commutative integrable equations derived from scalar Lax operators. 
The method can easily adapt to the case of (2 + l)-dimcnsional differential-difference integrable equation with 
scalar Lax operators |BS01j . 

Noncommutative integrable equations can be considered as quantised analogs of classical integrable equations, 
which is very important in string theory. Right after finishing the paper, we noticed a few papers on constructing 
noncommutative integrable equations in the framework of the Sato theory motivated by noncommutative gauge 
theories, cf. |HT03I IHam05| and refe rences therein. However, the authors only used an ansatz to derive the 
equations from some scalar Lax operators in case a instead of this systematical manner. There are few papers 
on noncommutative (2 + l)-dimensional integrable equations derived from the scalar Lax operators in cases 
b and c such as noncommutative mKP, where the new type of nonlocal terms such as (Dx + ^)~^Uy and 
{Dy — Cu)~^Ux appear. To study these equations further such as constructing their Hamiltonian operators and 
soliton solutions may provide a deeper understanding of noncommutative geometry. 
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Appendix: Proof of Proposition |H 

Proof. We need to show (cf. for alH > 2 that 

£=2 s=2 

^In the commutative case, this worked out in paper IBSOll by both central extension approach and the operand approach. 
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We prove it by induction. It is true for I — 2 since 



Assume that formula H28|) is valid for Z — 1. Then its the left-hand side equals 

L I c{is) x^l (is) Ai2) (12) A^i-2) Ah-2) -v- I <;(*i-i) I 

N, X^'^*") n(*=) n'^'^) . T7*^*'-i^ 



f(i2) (12) t(ii-l) I 



('^('2) „(i2) „(^i-i) TT^f^^'h 



h /ff'^) ^(^2) <;(il-2) „('!-2) Y , Ail--!-) , h ('f^*^) ^(13) V _J_ 



W ("^^f'^) (22) f(ii-l) „(«!-i) \- 1 i\r /fCil) Ail) Y^ 

/'t^*^) ^(i2) f(ji-2) „(«!-2) V u_ j_ 

iV;_2(?j, >^J^lUj, ,Aj 

^r Cf^^^) (12) A^i-i) Ail) 1 vaa^ rf^*') y\' 
Meanwhile, the right side gives us 

h ic(^2) (i2) Ail-2) Ail-2) -V- I I 

bi-sit'^t'' ■ ■ ■ ■ fc^ + ^gi^^^ + + 1') 



AT (C^^^^ Y\Ar (f^^^^ y^i 



The difference between the above two expression vanishes due to the induction assumption, which implies that 

AT /"f^'^) Ai2) Ail-2) Aii-2) -V- I f (*i-i)^f, (Ai2) (i2) Ail-2) Aii-2) v_Lf(*'-ih 

h rt^^s) ^(i2) ^(i,_2) I „(»i-2) I Y\ h (Ah) Ai'i) Ail-i) Y\ 



AT /^(ii-l) y\ l\r /'<r('2) (i2) Ail-l) Ail-i) Y\' 



By now we proved our proposition, o 
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